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We develop an analytical method for calculating local correlations in strongly interacting ID 
Bose gases, based on the exactly solvable Lieb-Liniger model. The results are obtained at zero and 
finite temperatures. They describe the interaction-induced reduction of local many-body correlation 
functions and can be used for achieving and identifying the strong-coupling Tonks- Girardeau regime 
in experiments with cold Bose gases in the ID regime. 



I. INTRODUCTION 



Recently, one-dimensional (ID) Bose gases have been created in long cylindrical traps by tightly confining the 
transverse motion of particles to zero-point oscillations P, S In the present stage, one of the main goals is to 
achieve the strong-coupling Tonks-Girardeau (TG) regime 0,|3|j where due to repulsion between particles the ID 
Bose gas acquires fermionic properties. These studies have revived an interest in dynamical and correlation properties 
of ID gases. In most cases, correlations and dynamics of trapped atomic gases in the ID regime can be investigated 
^ lL J.I A ^1 0> Q by using the Lieb-Liniger model ^ , which assumes that particles interact via a delta- function 
repulsive potential. Of particular importance are local 3-body and 2-body correlations [Til IT^ as they govern the 
rates of inelastic processes, such as 3-body recombination and photoassociation in pair collisions. The measurement 
of these rates provides a way to identify the strong-coupling Tonks-Girardeau regime, and the intrinsic decay process 
of recombination in 3-body collisions is crucial for the stability of the gas. 

The Lieb-Liniger model for a uniform system of N bosons belongs to the class of integrable models of statistical 
physics [T^ . The ground state and spectrum of elementary excitations for this model were found by Lieb and Liniger 
|5| and a theory for finding thermodynamic functions at finite temperatures was constructed by C.N. Yang and C.P. 
Yang [l^ . These quantities are determined by relatively simple integral equations for the distribution of quantum 
numbers (quasi-momenta or rapidities). On the other hand, the problem of correlation properties is far from being 
completely resolved, except for some correlation functions in the limiting cases of weak and strong interactions. For 
example, the case of infinitely strong interactions is to a certain extent equivalent to that of free fermions and the 
interactions play the role of Pauli principle 0|. In this limit, any correlation function of the density is given by the 
corresponding expression for fermions, and the one-particle density matrix is the determinant of a certain integral 
operator 15, 16] . The expressions for the one-body and two-body correlations for an arbitrary interaction strength 
were obtained by using the Inverse Scattering Method However, closed ana lyti cal results can be found only as 
perturbative expansions in the limiting cases of weak and strong interactions [l8l Il9l |2(]| . 

In this paper, we calculate local many-body correlations, that is the expectation value of a product of 2m field 
operators (to = 2, 3, . . .) at zero space and time separation. The locality allows us to obtain closed analytical results at 
zero and finite temperatures in the limiting cases of strong and weak interactions. The paper is organized as follows. 
In section m we present the Hamiltonian and introduce relevant parameters of the problem. Section IITTI contains the 
main result of the paper, the derivation of local correlation functions at zero and finite temperature in the strong- 
coupling regime. In Section Hvl we use this result for the particular case of two-body and three-body correlations. For 
the sake of completeness, local correlations for the weak-coupling regime, following from the Bogoliubov approach, 
are presented in Section We conclude in Section IVTl 



II. INTERACTING BOSE GAS IN ID 

We consider the system of A'' bosons on a ID ring of length L in the thermodynamic limit N, L ^ oo, with a fixed 
density n — N/L. The particles interact via a repulsive delta- function interaction potential and are described by the 
Hamiltonian 
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Here M is the particle mass, Xj is the coordinate of the j-th particle, and c = mg/h? ^ with g > being the coupling 
constant. The quantity c is the inverse interaction length for the 2-body problem with the delta- function potential 
g5{x). In other words, the wave function of two particles decreases on a distance scale 1/c as they approach each 
other. 

The Hamiltonian is diagonalized by means of the Bethe Ansatz 5] . The many-body wave function is symmetric 
with respect to permutation of particle coordinates, and in the domain < xi < . . . < xn < L the wave function 
reads 

<I>(xi , X2 , . . . , xat) = -1- 5^ a{P) S . (2) 
TV (c) ^ 

In Eq. Q the sum runs over A''! permutations P acting on N indices of quantum numbers (quasi- momenta) kj. These 
quantum numbers are chosen such that ki < k2 < ■ ■ ■ < kN ■ They determine the amplitudes 



of each term in Eq. |(2Jl and their quantization in a finite system follows from the requirement of periodicity of the 
wave function jSJ. The corresponding eigenenergy is E = ^k^ , and M{c) in Eq. 101 is the normalization constant. 
The latter has been calculated in Refs. pal2^. 

The key parameter of the system is the ratio of the mean interparticle separation 1 /n to the interaction length 1 /c: 

7 ==- = -2-- (4) 
n h n 

At sufficiently low temperatures, for 7 S> 1 the gas is in the strong-coupling regime: the many-body wave function 
strongly decreases at interparticle distances much smaller than 1/n. In the extreme case 7 = 00 the amplitudes (O 
are determined by the sign of the permutation: 

a{P) = i-lf. (5) 

So, the wave function Q vanishes when at least two particle coordinates coincide. In this limit any correlation 
function of the density operators can be calculated by using the theory of free fermions, since in the domain < 
xi < X2 < ■ ■ ■ < Xn < L the wave function ^ is the Slater determinant constructed out of the plane waves with 
momenta kj. In the case of large but finite 7, the calculations can be performed perturbatively as expounded in the 
next Section for the local density correlations. 

In the opposite limit, 7^1, the gas is in the weak-coupling regime. In this case the ground state of the system 
is well described by the Gross-Pitaevskii mean-field theory, and correlation functions at T = can be found on the 
basis of the Bogoliubov-Popov approach (see i21j]). This approach also covers a significant part of the weak-coupling 
regime at finite temperatures. 

For finite temperatures, aside from 1/c and 1/n, one has another length scale, the thermal de Broglie wavelength A = 
(27r?i^/mT)^/^. The latter is especially important at temperatures T > T^, where = h^n^/2m is the temperature 
of quantum degeneracy. The relation between the three length scales determines various finite-temperature regimes 
of the ID Bose gas, which we discuss below in the limits of strong and weak interactions. 

III. LOCAL CORRELATIONS IN THE STRONG-COUPLING LIMIT 

We now consider the strong coupling limit 7^1 and calculate the local m-particle correlation function 

g„(7) = ((*t(o))™($(0))™), (6) 

where '^'^ {x) (x) are the creation and annihilation operators of bosons. At T = the expectation value (. . .) is 
taken with respect to the ground state of the system, and at finite temperatures we assume the average in the grand- 
canonical ensemble. We first present the derivation at zero temperature and then generalize the method to the case 
of finite temperatures. 

For T = 0, in first quantization the correlation fvmction g,„ 10 reads 

gm(7) = UAT ' 77 d,Xrn+l---dXN % {0,...,0,Xm+l,---,XN) , (7) 
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where the ground state wave function ^q''''' is given by by the Bethe Ansatz expression (01. In the strong-couphng 
hmit the amphtudes a(P) |(2Jl can be expanded in inverse powers of the interaction strength 



<p) - n 



IC + Kp. 



ic — kp. + k} 



IC 



(8) 



The quasi-momenta are close to their ground state values at 7 = oo, given by the uniform Fermi-Dirac distribution 
in the interval [—kp, kp], with the Fermi momentum kp = irn. One then sees that in fact Eq. ((SJ gives the expansion 
in inverse powers of 7. For 7 3> 1, we extract the leading behavior of the wave function $q at m coinciding points by 
appropriately symmetrizing the amplitudes a(P) for each permutation: 



1 

— ^ a(Pp, ,Pp„...,Pp^, P^+i,. ..Pn)^ 7-w;;i3T)7iA™(/cp, ,...,kpj. 



^^^■Jm(m — 1)/2 

f 

where the sum runs over m! permutations p of numbers 1, 2, . . . , m, and 



(9) 



(10) 



i<j 



is the Vandermonde determinant. Up to the leading term in 1/c, the ground state wave function at m coinciding 
points is therefore given by 



<P^^\0,...,0,Xm+l,---,XN) 



(11) 



This allows us to express ^l^\o, . . . ,0,Xm+i, ■ ■ ■ ,xn) through spatial derivatives of the wave function of non- 
interacting fermions, 



<^^q°°\xu...,Xn) 



From Eq. Hll|l we find 

^^^\0,...,0,Xm+l,---,XN) 



m(m-l)/2 



Am{dxi,- ■ . ,dx^)^Q°°\xi, . . .,Xm,Xm+l, ■ ■ ■ , X n) 



(12) 



xi — ...—Xm—0 



Returning to Eq. Q , one sees that the local correlation function is given by derivatives of the m-particle correlation 
function of free fermions at Xi = . . . = x^ = yi = . . . = ym = 0- 



1 



(5^1, . . . ,dx^)/:^^{dy^,. . . ,dy^){il}\xi) . . .tpHxm)ip{ym) ■ . . V(2/i))- 



(13) 



The operators ■0^, ip are now fermionic field operators with canonical anticommutation relations. Using Wick's 
theorem we express the expectation value of these operator as a product of one-particle Green functions: 



(V'^(xi) . . .ip''{x,n)il^{ym) ■ --iPiyi)) = G{xi,yp^)G{x2,yp2) ■ . . G(x,m, ?/p„J. 



(14) 



The one-particle Green function is given by 

G{x,y)^G{x~y) = {^\x)^/j{y)) = i ^ 7V"(A;) e^'^^--^), 



(15) 



where iV°(fc) is the ground state Fermi-Dirac distribution. Then, from Eqs. H15|l . (|14|l and we obtain 



P 



ki ...k„ 



Xl 

yi 



(16) 
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where ^rn{dx) stands for Am(9a;j, . . . ,(?a:,„)- This expression can be further simpHfied as the Vandermonde deter- 
minant is a totally antisymmetric function of its variables. Therefore, reordering yp. arguments cancels the factor 
(— Using this fact, acting with derivatives on the corresponding arguments in the exponents, and noting that all 
permutations give the same result, we arrive at the equation: 

T)l4r E N\k,)...N^{h^)/^l{k^,...,k^). (17) 



(fn(m-l) ]^ 



The developed method was briefly outlined in Ref. mI'I for 3-particle local correlations. It is readily generalized to 
the case of finite temperatures. In the grand canonical ensemble the expression for the m-particle local correlation 
function g,n © reads: 

^ lTr{(vI/t(0))"(vI/(0))'"exp{-(ff-/i7V)/r}} = 

N=l 'a 

where Z is the partition function, z = eKp(fj,/T) is the fugacity, ^ is the chemical potential, and the index a labels 

eigenstates of the system. Expanding the eigenfunctions <I'L'^'' at coinciding points exactly as above in the case of the 
ground state wave function, we express the leading contribution to g„ through derivatives of the finite-temperature 
correlation function of free fermions ai xi — X2 — ■ ■ ■ ^ x„i = yi = . . . = ym = 0: 



9n 



'^midx) \n{dy) {ijj'' i^l) ■ ■ ■ 1pHxm)lJKym) ■ . . ■(/'(2/i))t (19) 



The finite-temperature m-fermion correlation function is again calculated with the use of Wick's theorem. This gives 
Eq. (|14|l in which G{xi, yp.) are now finite-temperature Green functions for free fermions. We then proceed in the same 
way as at T = and arrive at Eq. H17|l . with N'^{k) replaced by the Fermi-Dirac distribution at finite temperatures: 

AKL.^ ^exp (-A2fcV47r) 

"^(^^ = l + zexp(-A^fcV4^) - (2°) 

At temperatures T ^ the characteristic momentum of particles is the thermal momentum fc^ ^ 1/A. Therefore, 
the small parameter for the expansion of the amplitudes a{P) in Eq. (jSl is 1/Ac. Thus, one must satisfy the inequality 
Ac ^ 1, which requires temperatures 

T < -f^Td. (21) 

From this point on we consider together the cases of zero and finite temperature. In the thermodynamic limit, 
expressing the momenta in terms of the size of the Fermi zone, k = 2kFX = 2'Knx, we have 

27r \ ^'^(^^1) 

m! (An, z) , (22) 

where the m-fold integral /,„(An, z) depends on the ratio of the de Broglie wavelength to the mean interparticle 
separation and on the fugacity z (or chemical potential), and is given by 

/ + 00 
dxi .. . dxmN{xi) . . . 7V(xm)A^„(a;i, . . . ,Xm), (23) 
-oo 

where N{x) is given by Eq. (|20|l with k = 2'Knx. In order to keep the density constant one must fix z as a function 
of An by the normalization condition 

-t-oo 

N{x)dx = 1. (24) 

-oo 

Under this condition the integral in Eq. (|23|1 becomes a function of An only. For calculating this function we employ the 
method of orthogonal polynomials used in Random Matrix theory [23 . Consider a set of polynomials Pj {x) — x^ + . . . 
for j = 0, 1, 2, . . ., orthogonal with the weight N{x) on an infinite interval: 

N{x)P^{x)PJ{x)dx = h^S^j. (25) 



Then the value of the integral (|23|l is expressed through the normalization coefficients hj as follows: 

Im {■^n,z) = mlho . . . hm-i- 



(26) 



At T = we have An = oo. In this limit the distribution function N{x) is uniform in the interval — 1/2 < x < 1/2, 
and is zero otherwise. Polynomials -Pj(x) can be expressed in terms of Jacobi polynomials and we obtain 



/,„(An) = m! J| 



[r(j + i)r 



r(2j + i)r(2j + 2)' 



An = oo. 



From Eq. H22|l we then obtain the local correlation function 



.(7) 



m(m — 1) 



where the coefficients Am satisfy the recurrence relation 



A. 



m+l 



2m + 1 

22m+2 



r(m + 2) 



r(rn + 3/2) 



Ai = l. 



(27) 



(28) 



(29) 



Although general expressions for the normalization coefficients hj in Eq. (|25|l are not known, for sufficiently small m 
these coefficients can be calculated straightforwardly as they are related to the moments of the distribution function 
N{x). For example, one finds 



/lo = 1, hi =< >, h2 =< x'^ > - < x^ < x^^ >= 



x^^N{x)dx. 



(30) 



We finally consider high temperatures T ^ Td that still satisfy Eq. (|21|) . In this case the fugacity is 2; = An, and 
the distribution function is Gaussian: N{x) — z exp (— TrA^n^a;^). The polynomials Pj{x) are Hermite polynomials, 
and we find 



ml 



I„i(An) = — , , , , , 

^ ' (\/2^An)'"(™-i) fj-^ 

Then Eq. (|22ll for the local correlation function takes the form 



iynr(l+i), l/c«A«l/r 



9n 

n" 



Bn 



m(m — 1) 



Ac 



1/c < A < 1/n, 



(31) 



(32) 



where the recurrence relation for the coefficients Bm reads Bm+i = (m + l)r(m + 2) B„n and Bi = 1. 



IV. THREE-BODY AND TWO-BODY CORRELATIONS 



In this Section we use the general results of Section IIIII for the particular case of three-body local correlations 
in the strong-coupling limit. Two-body correlations have been discussed in Ref. [T^ and are expounded here for 
completeness. 

At temperatures T ^ we have An ^ 1, and the local correlation functions 32 and 173 are close to their zero 
temperature values. Thus, the system remains in the Tonks-Girardeau regime. For calculating 52 and we use 
Eqs. H22I) and H2t)|l . with normalization coefficients hi and /i2 following from Eq. H30|) . The quantities < > and 
< > at finite T are obtained on the basis of the Sommerfeld expansion: 

2, ^ ^ 2, . 2 2 j 1 I 2j (T " 



< >. - < X- >o^ = y- \ . (33) 



This gives the following expressions: 
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The results of Eqs. (|34() and (|35|l have a clear physical explanation. For T <^ Td the characteristic momentum of 
particles is of the order of kp = T^n. Fermionic correlations are present at interparticle distances x ^ 1/c, since for 
smaller distances the correlations do not change. Therefore, 32 is nothing else than the pair correlation function for 
free fermions at a distance x ^ 1/c T^]. This function is of the order of [kp /cY ~ {tt/"/)"^, which agrees with Eq. 134|l . 
With the same arguments, one finds that the 3-body correlation function is 173 ^ (kF/cY ~ (7r/7)^, which coincides 
with the result of Eq. Ip!5|l . 

For temperatures in the interval ^ T ^ 7^Tk, where 1/c <C A <C we use Eq. (|32|l directly and obtain: 

(36) 



52 


Stt 


2 T 




(Ac) 2 






5767r3 


9 /T 




(Ac)6 





(37) 

This regime can be called the regime of "high-temperature fermionization" The gas is no longer degenerate, 

but due to strong interaction between particles the correlations still have a fermionic character. With regard to local 
correlations, the main difference from the low temperature case is related to the value of the characteristic momentum 
of particles. It is now of the order of 1/A, instead of fc^ at T ^ T^. Accordingly, the two-body and three-body local 
correlation functions are g2 ~ (l/Ac)^ and 53 ~ (1/Ac)^, as described by Eqs. and l|T7|) . 



V. WEAK-COUPLING LIMIT 



In the weak-coupling limit, where 7 <C 1, one can rely on the mean- field approach. The mean-field interaction 
energy per particle is proportional to ng and it is reasonable to introduce the correlation length Ic = h/y/mng. At 
temperatures T <^ T^, over a wide range of parameters the correlation length 1^ l^, where is the phase coherence 
length. Then the equilibrium state is a quasicondensate, that is the state in which density fluctuations are suppressed 
but the phase still fluctuates [2^. A review of earlier studies of ID Bose gases in the weak-coupling limit may be 
found in Ref. [2l] |. The zero temperature phase coherence length is always larger than Ig., and at finite T one has 



h n/mT. Therefore, the condition Ic ^ can be rewritten as 



(38) 



One then sees that the quasicondensate regime requires sufficiently low temperatures or relatively large interaction 
between particles. 

Due to the presence of phase coherence on a long distance scale, local correlation functions in the quasicondensate 
regime are close to n^. Deviations from this value can be calculated straightforwardly by using the Bogoliubov 
approximation [27j . For 172 the results are known (see [T^ l28j and references therein), and 33 at T = is given in 
Ref. [ijLHere we complete the picture and present general results for gm ■ They are obtained in the same way as 32 
in Ref~[l3. 

We represent the bosonic field operator as a sum of a macroscopic component ^Pq and a small component ^f' 
describing finite-momentum excitations. To be more precise, the component ^'o contains the contribution of excitations 
with momenta k ^ kg <^ whereas 5*' includes the contribution of larger k. At the same time, the momentum ko 
is chosen such that most of the particles are contained in the part 'Sq. This picture is along the lines of Ref. [2ll |. 
and the momentum ko drops out of the answer as the main contribution of the excitation part \E'' to local correlation 
functions is provided by excitations with k Confining ourselves to the terms quadratic in and taking then 



into account that l^'o 



2m 



^i(\I>'t\j/'^^ the m-particle local correlation function is given by 



m(m — 1) / , ,+ „ 



(39) 



The normal and anomalous averages, (^P'tv]/') and (^'^'), can be calculated by using the same Bogoliubov trans- 
formation for as in the 3D case: 



,Tf/ \ ^ / ikx t —ikx 



(40) 



where ak and aj. are annihilation and creation operators of elementary excitations, Uk,Vk — {sk ± Ek)/2\/ekEk, 
£k = -vZ-El + 2ngEk is the Bogoliubov excitation energy, and = fi'k^ I2m. This immediately gives 



5m , , , dk 

= i + mm— 1 / 

^ 27m 



El. 

£k 



(41) 
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with Nk = [exp (sk/T) — being the equilibrium occupation numbers for the Bogohubov excitations. 

The integral term in Eq. H41I) contains the contribution of both vacuum and thermal fluctuations. The former is 
determined by excitations with k ^ and at T = we find 

Qm m(m — 1) _ , , 

^ = 1 - ^ -Vl- 42 

At temperatures T <^ ng ^ jTd, thermal fluctuations provide only a small correction {7ry/^/24:){T/jTd)'^ to the result 
of Eq. 

For temperatures T 3> ng ^ ^Td (but still T <C ^J^Td), the contribution of thermal fluctuations to the integral term 
in Eq. H41|l is the most important. It comes from excitation energies Su ^ ng and, hence, one may put N{k) — T/sk 
in the integrand. This yields 

Qm m(m — 1) T , ^ 

— = 1 + — ^ '- . (43) 

4 ^Td ^ ' 

One clearly sees from Eq. (143(1 that gm increases with decreasing interaction strength or increasing temperature. 
The ratio T/Td^J^ is simply lc/l,j,- For Ic ^ or T ^ Tdy^, the gas leaves the quasicondensate regime and Eq. 
is no longer valid. At a temperature and interaction strength satisfying the condition T ^ Tdy/^ the gas is in 
the decoherent regime where the interaction between particles is much less important. Note that for a very small 
interaction strength this regime is present even at ternperatures far below Td- The cross-over from the quasicondensate 
to the decoherent regime was discussed in Refs. 0,113 on the basis of calculations for g2- In the decoherent regime 
the local correlation function is close to the result for an ideal Bose gas following from Wick's theorem: (?m/ri™ = ml. 
The corrections to this value can be calculated perturbatively 12]. The lower is T, the smaller is 7 at which the gas 
enters the decoherent regime. For T = this transition is discontinuous and occurs at zero interaction strength. 



VI. CONCLUDING REMARKS 



From a general point of view, local correlations of the ID system are less universal than the long-wave behavior. 
Local correlation functions depend on a particular model used for calculations. Our approach leads to physically 
transparent analytical results for g^i of the ID Bose gas in the Lieb-Liniger model. 

The use of this model for trapped atomic gases in the ID regime is justified by the short-range character of 
interatomic interaction: the characteristic radius of interactions is much smaller than any length-scale of the Lieb- 
Liniger model. For an atomic gas in an infinitely long cylindrical trap, harmonically confined with frequency luq in the 
transverse direction to zero-point oscillations, the coup ling constant g is expressed through the 3D scattering length 
a and reads '2^ g = ih^a/Ml"^, where Zo = Mwq is the amplitude of zero point oscillations. Then, to describe 
the system by the ID Hamiltonian (QJ, it is sufficient to satisfy the inequalities Iq « 1/ri, At [Tllll^. 

The reduction of 53 in the strong-coupling limit, which follows from our results, is important in two aspects. First, 
it indicates a possibility of achieving this limit at a high gas density (large number of particles), since the rate of 
decay due to 3-body recombination is proportional to 173 and will be strongly suppressed. Second, the measurement 
of 3-body losses of particles can be used for identifying the strong-coupling Tonks-Girardeau regime and the regime 
of high-temperature fermionization. The identification of these regimes can also be provided by the measurement of 
photoassociation in pair atomic collisions as the rate of this process is proportional to the two-particle local correlation 
function 52 • In this respect, it is worth mentioning that recent Monte Carlo calculations of 172 at T = for the number 
of particles as low as 100 |33 agree with our results obtained in the thermodynamic limit. 
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